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Generalized Lovelock gravity 

Qasem Exirifard 

School of Physics, Institute for Studies in Theoretical Physics and Mathematics (1PM), P.O.Box 19395-5531, Tehran, Irai^ 

In the Riemann geometry, the metric's equation of motion for an arbitrary Lagrangian is succinctly 
expressed in term of the first variation of the action with respect to the Riemann tensor if the 
Riemann tensor were independent of the metric. Let this variation be called the E-tensor. 

Noting that the E-tensor and equations of the motion for a general Lovelock gravity have the same 
differential degree, we define generalized Lovelock gravity as polynomial scalar densities constructed 
out from the Riemann tensor and its arbitrary covariant derivatives such that they lead to the same 
differential degree for the E-tensor and the metric's equation of motion. 

We consider Lagrangian densities which are functional of the metric and the first covariant deriva- 
tive of the Riemann tensor: £(V a Rijki,gij). We then present the first non-trivial examples of the 
generalized Lovelock gravity terms. 

Quantization of gravity requires the Einstein-Hilbcrt action to be corrected. But what are the form of the cor- 
rections? In string theory, corrections in the form of the covariant derivative of Riemann tensor exist in the sub- 
sub-leading levels [1-3]. In the Horava-Lifshitz proposal for quantum gravity in 3 + 1 dimensions Q, the quantum 
corrections are functional of the covariant derivative of the Riemann tensor. So quantum corrections in d = 4 most 
probably will depend also on the covariant derivatives of the Riemann tensor. 1 This dependency, however, increases 
the number of possibilities for the Lagrangian density: the number of algebraically independent polynomial scalar 
densities constructed out from the Riemann tensor and its covariant derivatives are larger than the number of those 
constructed out from the Riemann tensor alone. Thus it seems plausible to provide a criterion in order to systemati- 
cally and self-consistently selects a subset of the Lagrangians densities from all possibilities, a subset that also includes 
the covariant derivatives of the Riemann tensor. In this note, we aim to provide such a criterion. 

A natural rephrasing of this aim is what the generalizations of the Einstein-Hilbert equation are. David Lovelock 
long time ago posed this question. He demanded the generalized Einstein tensor to be: 

1. divergence free. 

2. symmetric. 

3. concomitant of the first and second derivative of the metric. 

In a set of publications (llT-fl3| , he then proved that the only Lagrangians meeting his demands are 
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1 In the phenomenological geometric approach, the dark energy and dark matter problems are evaluated as needs to modify the dynamics 
of the space-time in very low space-time curvature. The simplest geometric models to resolve these 0,01 EH rely in the f(R) gravity ||. 
However these models do not reproduce the Tully-Fisher relation [8| and they introduce a parameter that depends on the (baryonic) 
mass of the galaxy in order to account for its observed anomalous rotational velocity curve. Inclusion of the covariant derivatives of 
the Riemann tensor into the action resolves these problems [9( . So studying actions that depends on the covariant derivatives of the 
Riemann tensor is admirable even within the phenomenological geometric approach to dark paradigms. 
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where [-j] represents the integer part of -j, and a„'s are some constant values. Neither of the above Lagrangians 
includes the covariant derivative of the Riemann tensor. So we do need to rephrase and modify the Lovelock criteria 
such that further possibilities are allowed. 

Having noted that "a characteristic feature of Lovelock terms is that their first non vanishing term in the expansion 
of the metric around flat space-time is a total derivative, S. Cnockaert and M. Henneaux have investigated generalized 
Lovelock terms defined as polynomial scalar densities in the Riemann curvature tensor and its covariant derivatives 
(of arbitrarily high but finite order) such that their first non vanishing term in the expansion of the metric around 
flat space is a total derivative " [l4| . They however have reported that their generalized Lovelock terms contains only 
the usual ones. 

Narcsh Dadhich has provided a new independent identification of the Lovelock gravity from the Bianchi derivative of 
a curvature polynomial [TBj . One can investigate if this classification leads to a non-trivial generalization of Lovelock 
gravity. Ref. [lg, 113] investigate and report the consistency 2 of the Palatini (first order) (l8l - l20j formulation and the 
metric (second) formulation of the Lovelock terms. Ref.|l6|, having evaluated the consistency of the first and second 
order formulations as a criterion to restrict the form of the Lagrangian promises to apply this criterion on other 
Lagrangians. This promise is yet to be fulfilled. 

In this note we would like to look at a characteristic of the Lovelock gravity that so far has been overlooked in all 
attempts to generalize Lovelock gravity terms. To illustrate this characteristic let us look at the Gauss-Bonnet term 
which is n — 2 term in ([T]) 

L GB = - 4Rtf J# + R 2 , (5) 

The equation of motion for each of R^ kl fjym, RijR 1 ^ and R 2 is a fourth order differential equation. But the Guass- 
Bonnet term is the combination of these terms that brings down the degree of differential equation of motion by two. 
Let us examine if the degree of the equation of motion can be brought down by two for a general action. In so doing 
we need the equation of motion for a general action. 

A general action for the metric of a D dimensional space-time can be presented as follows 



S, 



gravity 



J d D x^/-dctgL[g ij ,R ijk i,V i } , (6) 



Llgi^RijM, Vi] = C[ 9ij , Ri jkl , V a i? yfc \ V (oi V a2) i? yfc( , • • • , V (ai • • • V a „)i^ fci ] , (7) 

where n is a natural number. The first variation of the action with respect to the metric then gives the equation of 
motion for the metric 

3C 1 

-T» = — + -Cg ij +E\ M W a ?~< + 2V a V p E ia ^ , (8) 

Em = d£__ Va d£_ +t+{ _ dc (9) 

OUijkl CV a-ttijkl " V (ai ' ' ' V a n )-tHjkl 

where T lJ is the energy momentum tensor of matter fields minimally coupled to gravity, and partial derivatives of L 
are taken assuming that gij, Rijki > • " > V( 0l • • • S7 an ^Rijki are independent variables, and partial derivative coefficients 
are uniquely fixed to have precisely the same tensor symmetries as the varied quantities [2l[ . The E-tensor 3 would 
be the first variation of the action with respect to the Riemann tensor if the Riemann tensor were independent of the 
metric. 

In the equation of motion ©, VaV^fi 1 "^ generally leads to the appearance of two extra derivatives of the metric 
which are absent in the first three terms of the r.h.s of ((SJ. In other words the differential degree of the equation 
of motion is generally two degrees higher than that of the E-tensor. Requiring the same differential degree for the 
E-tensor and the equations of motion, thus, is a criterion to single out a specific set of Lagrangian from a larger 
given set. We refer to this criterion as the E-criterion. 4 The chosen/restricted subset has the privilege of leading to 
a differential equation of a lower degree than that of a general Lagrangian. Lovelock gravity terms are examples of 



2 This is a strong form of the consistency, in the sense that all the solutions of the metric formulation are also the solutions of the first 
order (Palatini) formulation. Ref. 1 1 6t) shows these two formulations are equivalent in asymptotically 'flat' space-time geometries. 

3 Ref. [2^ | uses P to refer to this tensor. We use the notation of [2l| and call it the E-tensor. 

4 In order to have the Einstein-Hilbert action included, we should have rephrased the E-criterion as terms for which the E-tensor's degree 
is not larger than the degree of the equations of motion . We are, however, interested in the higher derivative corrections/modifications. 
So requiring the same degree for the equations of motion and the E-tensor suffices for our purpose. 
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these terms. So we name terms that satisfy the E-criterion as the generalized Lovelock gravity terms. In the next 
section we consider Lagrangians which are functional of only the first covariant derivative of the Riemann tensor. We 
obtain a sufficient condition for the E-criterion. We find a family of solutions for this condition. We then discuss on 
the uniqueness of the these terms. At the end, we will provide the summary and outlooks. 

I. FUNCTIONAL OF THE FIRST COVARIANT DERIVATIVE OF THE RIEMANN TENSOR: 

A Lagrangian in the form of C[gij, V a Rijki], polynomial in terms of V ' a Rijki is a summation of the following terms 

where C ai "' 0nll "' i2 '' jl "'- ?2n is a functional of the metric components, and n is a natural number. Since 
£oi~ a„n~ %2„,h---hn j s a functional of only the metric's components it holds 

VpC oi-o„ii-i 2 „ii-i2„ _ o. (11) 

and it carries an even number of indices. Therefore n is a natural even number. 

We note that only the part of C 0l '" 0n * 1 '" t2nJ ' 1 '"-' 2n which is symmetric under the exchange of 
(a m ,i2m.-l,i2m,hm-i,hm) with (a p , i2p-i, hp, hp-x, hp) for all m and p contributes to the Lagrangian density. 
So we choose C^'" " 11 '"* 2 ™- 31 '"^ 2 ™ such that it satisfies 

Qa 1 ---a p ---a m ---a n i 1 ---i 2p - 1 i 2p ---i2rn.-li2rn---i2njl---j2p-lj2p---j 2 m-lj 2 m---j2n 

Let it be recalled that all the components of the Riemann tensor and its covariant derivatives are not algebraically 
independent. The Riemann tensor constructed out from the Levi Cevita connection satisfies: 

(13) 
(14) 
(15) 
(16) 

Because the Lagrangian density is the multiplication of the Riemann tensor and the C-tensor, it can be deduced 
from (TiUj) , (fLl|) , and (Ti"5)) that only the part of the C-tensor that owns the following properties contributes to the 
Lagrangian: 







= o, 






= o, 


Rili23lj2 ~ 




= o, 


■i\i2]2a\ + Via 


^Hliaaiii 


= 



Qa 1 ---a n i 2 ii---i2 n Ji---]2n _ " '%» l»2 • • -12n3l • • • J2n (17a) 

Qa 1 ---a n i 1 ---i2 n j 2 ji---j2n _ _Qa 1 ---a n i 1 ---i 2n jij 2 ---j2n (17b) 



^ai'"a„Jlj2»3— »2n»l»2j3--j2n _ ^7ai"-O n »l — I2n31 — J2n (17c) 

The above relations are nothing more than saying that when an scalar is constructed out from the direct multiplication 
of a symmetric tensor and a general tensor, then only the symmetric part of the general tensor contributes to the 
scalar. We choose C such that it explicitly holds (jl7l) . 
Now let it be noticed that (Til?]) implies 



/70i— o n *l— ianJl— janfy /,> . V /.' _|_ V 7? ; (| ' I .O, 

rewriting which yields 

^Oi02— <»nili2i3—»anil— J3n _|_ Qha2—a n i2aii3—i2njl—j2n _|_ ^202— CSraCtljlj3— i2tiil-"i2n ^y R =0 (19) 

which implies that only the part of the C-tensor that holds 

Qaia 2 ---a n i 1 i2i3---i 2 nji---]2n _|_ Qi\a,2 -- a n i2aii3 — i2nj\— 32n _j_ £»*2<»2 — a n aiii 13 — i2njl— J2tl _ Q (20) 
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contributes to (|10[) . In order to further clarify this statement, let the C-tensor be re-expressed as follows 

ga 1 a 2 ---a n i 1 i 2 i3---i 2n j 1 ---j 2n _ ^ja 1 a 2 ---a n i 1 i 2 i 3 ---i 2n j 1 ---j 2n _j_ ^a 1 a 2 ---a n i 1 i 2 i3---i2njl ■■■jin (21) 
3y^ a i°2 ■■■a„iii 2 i3 ■■■i 2n jl ■■■j 2 n _ Qaia 2 ---a„iii 2 i3---i 2n ji---j 2n _|_ £>ii 02 ■ ■ -a n i 2 a\i3 ■ ■ -i 2 „ ji ■ ■ ■ j 2 „ _|_ Qi 2 a 2 ■■ -a„aiii is ■■ -i 2n ji ■■ -j 2 „ (22) 



Then when we insert (f2lj) into the Lagrangian density, due to (|T9|) . we see that only the (7 a i a 2 - a n 4 i 4 243 -'2nji - j2n _ 
the part of the C-tensor that holds (|20|) - contributes to the Lagrangian density. So we choose the C-tensor such that 
it meets ([20]) . Note that neither this choice nor (fT2"|) . nor (fTT)) affects the generality of the considered Lagrangians. 

Eq. (TIT)) and indicates that the C-tensor carries all the symmetries of the covariant derivative of the Riemann 
tensor. So the partial derivative of the C n with respect to S7 ai Rijki can be simply written by 

IT n 

" _ r ,f>a.ia 2 ---a n i 1 i 2 i 3 ---i 2n ji---j 2n T~V i-j p. ("oo - ! 

Qry p — " U 11 V a p rL l 2 p-ll 2 pj2p-lJ 2 p ■ 

(JVa-tli 1 i 2 j 1 j 2 p=2 

If we had not chosen the C-tensor to satisfy (jl2"1) . (IT71) and (f^Uj) . then 3V could not been written in the above 

compact form. The E-tensor can be written in a compact form too: 

- C/ £„ 1 r)T7 R _ ^ ' 1 a 2 rl »3»4j3j4 J_J_ V a p - n 'J2p-lJ2pj2p-lj2p • 

»VoftiiI 2 jlJ2 p=3 

Note that the E-tensor is a functional of the Riemann tensor and its first two covariant derivatives. The equation of 
motion for £ n also can be simplified to 

r) 1 \ ™ 

+ ) JI V apJ R l2p _ ll2p , 2p _ lj2p + ET^R\ M + 2V a V p E a c f b = , (25) 

for which the E-criterion requires V Q V ' pE a ^^ b to be a functional of the Riemann tensor and its first two derivatives. 
From (HQ), V^Vj^ 1 "^ 2 follows 

V 'J2 V Jl- C/ £„ — V »2 V Jl V a! ™ p l^OJ 

f v a I >-i 1 i 2 j 1 j 2 



which can be re-expressed by 

^i^hE l C nn = -(ViJV^.V,,,] + [V^VaJV^+V^ViJ— — -2 , (27) 



Noting that commutators can be expressed in terms of the Riemann tensor, we conclude V; 2 V jl E% afih and 

V a , Vj, Vi 2 av7 p £ " have the same differential degree. So applying the E-criterion on C n is the same as requiring 

^o, V,-, Vi, flV7 p £ " to be a functional of the Riemann tensor and its first two covariant derivatives. 



Now let us look at V,„ V 



'2 v ji av ai fl, 1I2J1J2 - 

V- V- ^aiQ2---a„iiJ 2 43---i2„Jl---j2„Y7. v I TTv i? I (281 

2 Jl ^)V7 R _ V »2 V Jl I J_J_ V Qp- fl 42p-lJ2pj2p-U2p I ! k ZO J 

a V a-fl.li i 2 JU2 



\p=2 / 

expanding which yields 

V- V- - — n(n - 1 )C a i a 2---o. n iii 2 i3--->'2 n ji---j2 n y. V - V 7? TT V /?■ ■ ■ + 

2 31 r)X7 R ~ ^ ' 22 Jl a 2-''' l 3Uj3j4 J_ J_ V a p I -M2p-ll2p]2p-l]2p ' 

OV a-tli 1 i 2 j 1 j 2 

+n(n — l)(n — 2)C ai ° 2 a ™ IlJ 2*3 » 2 „ji j 2 „ \7 a2 Ri 3 i 4 j 3 j 4 Vi 2 V a3 i?i 5 i 6: , 5: , 6 x 

V ' apRi 2 p-ii2pj2p-lj2p ■ (29) 

p=4 
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"fan 


[^i 2 j^a 2 ]Ri 5 i e j 5 j e j 


Qa\a^ ■ ■ -a n 


ili2»3-"»2»j 


.31" 


■fan 


v 42 


v* 


V Q2 




_ ^ai«2- 


•a? 


1*1 *2 »3 - 


••»a> 




"fan 


[Via; r^jii ^a 2 ]] Ri 3 iij 3 ji 



Now suppose that it holds 

/»ai02"-On<ii2i3—»2n3i— 22n Y7. R. fl /Qm 

^ V I p -n.l2p-lt2p32p-l32p — U ; V JU y 

then 5 

(31a) 
(31b) 
(31c) 

Inserting (J3TJ) into (J29J) yields 

«/• " 

V- V- — - — n ( n _l)raia2"-a n iii2i3"-i2nji-"hn. (\^T7. [y . y 11 d. . . )TTn n . , 

2 J1 flV D ~ ' V" J J L ' U V *2! L V J1J V a 2 \\ rl t3U333i) 11 V a p /1 '*2p-l*2pj2p-lj2p ^ 

t/V a-Ki 1 i 2 j 1 j 2 p=3 

+n (n-l)(n-2)C aio -" a » w ^-W 1 -i2 n ([V 3l , V^]^^)) ([V, 2 , V^]^^) x 

n 

^ apRi 2 p-ii2p32p-lj2p ■ (32) 

(f3"2")) proves that when (|50|) holds then Vi 2 Vj x dv ^ Cn is a functional of the Riemann tensor and its first covariant 
derivative. So when (|30| holds V ttl V^ 2 Vj x 9V is a functional of the Riemann tensor and its first two covariant 

derivatives: the E-criterion for C n is satisfied. Recalling ^[i 1 Ri 2 i 3 ]faj 3 = 0, (|30|) is guaranteed to be satisfied if 

Qa 1 ---a n iii 2 i 3 ---i 2n j 1 ---j2n _ Qai ■ ■ -a„i 2 i 3 ii • • -i 2 „ fa ■ ■ -j' 2n _ ^jai •••a„i 3 iii 2 ■ ■ ■i 2 njl —fan (33) 

Let it be highlighted that ([33| is a sufficient condition for the C-tensor to meet the E-criterion. 



A. Finding the first generalized Lovelock gravity terms 

The 9 tensor (@J yields 



V ip J, -»2p-l»2pj2p-lj2p ~~ U ; V°V 



2p- £L 22p-l22pj2p-lj2p 
aili2l3—»2nj'l—j2n — fltaisil— fenil— js» fl ! 3 il h ' ' 'ian Jl ' ' 'J2 



(35) 



which forgetting the ai---a n are similar to (1301) and ([33]). This similarity motivates us to consider the following 
Lagrangian as an ansatz that meets the E-criterion: 

r _ (•ai—a n mi—i2 n ji—32n y /> y p. . . /"3fi\ 

^ u v ai - tL «l«2jlj2 V "i. iL *2?l— l*2n.72»-l.?2n ' V OU 7 

where 6* is given in (|4]) and C ai "' a ' 1 can be chosen to be symmetric under exchange of each of its two indices. 
We note that though e«i"-<*»0H-"Wi--j2n yields $8$ and (|33j) : 

(jai— a n Qhhh— i2njl"-fan — Qa,i—a n Qi2i3H-"i2n3l— fan — (jai— an ^31112 —i2njl-"fan (37) 



it does not have all the symmetries of the C-tensor: it does not hold (|20|) . We should pay attention that (|33[) is a 
sufficient condition for the E-criterion provided that the coefficients of the covariant derivative of the Riemann tensor 
hold all the required symmetries of the C-tensor. 

Let us decompose {:•'■—'■ I) ••• h-fa, to tne part tnat mee ts (JT2J), (JT7]| and (HOI), and the rest: 

Qa X "-a n Qii---i2n3l"-fan — £K»l" - a»il" - *3nil" - .j3n _(_ _ ( 4<»l - "<»n*l" , *3nJl"\j2n ^3g\ 



5 Commutators are used to have Vi 2 and act on the Riemann tensor before V a . 
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where C ai '" anil '" i3njl "' jan respects ([12]). ([TTf and (|20]). while ji 4. ai "" a " <1 "" ia, » J ' I '" J ' an does not respect these symmetries. 
It is worth noting that only the C' ai "' a ™ ll "'* 2 '»^ 1 "' J2n contributes to the Lagrangian density 

f — /'ai— Onflil— *2njl—jan Y7 E> Y7 U ( Qq\ 

/*(01— Onil — »2nil— J2n V /' V7 /> . . 

^ v ai Jt 2l22jlj2 V a n ±L l2n-ll2n32n-l32n 7 

Lagrangian density is independent of A 01 "'"" 11 "' 12 "^ 1 "'^ 2 ". It then follows from (l37t and (|38|) that the 
Qa 1 ---a n i 1 ---t2njx---02n meets ([55]). 6 So ([M]) meets the E-criterion and indeed is a generalized Lovelock Gravity. 



II. ON THE UNIQUENESS OF THE FOUND GENERALIZED LOVELOCK GRAVITY TERMS 

In this section we wish to test/prove uniqueness of some special cases of (|36|) . Let us consider the simplest choice 
of C n , the first example of for n = 2: 

C 2 = V a Ri jk iV a R ijkl - 4V J2y V°i^ + V a RV a R . (40) 



which should be called the generalization of the Gauss-Bonnet Lagrangian density. 7 ' 8 Perhaps (l40l) is not the most 
general Lagrangian quadratic in term of the first covariant derivative of the Riemann tensor. But is it the only one 
satisfying the E-criterion? 

In the following, we shall prove that f|40[) is the only Lagrangian quadratic in term of the first covariant derivative 
of the Riemann tensor meeting the E-criterion. In so doing we notice that the algebraically independent scalars which 
are quadratic in the Riemann tensor and quadratic in the covariant derivative read [25j 



RDR, V p V q RR pq ,V s VrR Pq R prqs , R pq UR pq (41) 

W P RV p R, V r R Pq VrRpq,VrR Pq V g Rpr ) V* R PqrS V tRpqrs ■ 



The most general Lagrangian density which is quadratic both in the Riemann tensor and the covariant derivative is 
a linear combination of all the above possibilities 

L = Cl RDR + c 2 V p V q RR pq + c 3 V r V s R pq R prqs + c 4 R pq aR pq + (42) 
+c 5 \7 p RV p R + c 6 V r R pq V r R pq + c 7 \7 r R pq \7 q R pr + R pqrs V t R pqrs , 

where c\ ■ ■ ■ c% are some constants real values. Performing integration by parts and using 

V p i? P9 = \^ q R, (43) 
V s R prqs = V p R qr - V r R qp , (44) 
the Lagrangian density (|42p can be rewritten to 



L = (c 5 - y - Cl )V p RV p R + (c 6 - c 4 - c 3 )V r R pq V r R pq - (c 7 - c 3 )R pq V r V p R qr + c^tR^V'R^ . (45) 



Noticing the algebraic identity of 



V r V p R qr = [V r ,\7 p }R qr + V p V r R qr = [V r ,\7 p }R qr + ^V p V q R, (46) 



e Note that also A a i-"" Xnil -'" i ^n-33 n has the symmetries of p3) . 

7 In contrary to the Gauss-Bonnet Lagrangian 1 140 I I is not a topological term in D = 4. This can easily be verified by evaluating the 
action of C2 for : 

dr 2 

ds 2 = -A(r)dt 2 H hr 2 (de 2 + sin 2 8d4> 2 ). 

A(r) 

The £2 action evaluated for J7{ reads: 

S C2 = 4tt J dt J dr^-((-4A' + 4A"r - 2A"'r 2 )(A - 1) + r a A" 2 + 3A' 2 + r 3 A"'A' - 4r 2 A"A') , 

which is not a total derivative or a topological term: ^ 2 5^ 0. 

8 The explicit form for n = 4 can be written using the compact form of the forth order Lovelock gravity that is presented in [23| . 
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and performing an integration by part, and using (|43p . (I45[) can be re-expressed by 

L = (c5-y- Cl + ^^)V p i?V p i?+(c 6 -C4-C3)V,.i? P9 V^^ (47) 

The commutators of the covariant derivatives can be expressed in term of the Riemann tensor. So R pq [\7 r , V p ]R qr 
is cubic in term of the Riemann tensor. 9 We are interested in actions which are quadratic in term of the covariant 
derivative, we thus set cj — c 3 . We also redefine the constants values 

c 2 

ax = c 5 - — - ci , (48) 

a 2 = c 6 - c 4 - c 3 , (49) 
a 3 = c 8 . (50) 

The most general Lagrangian density which is quadratic both in the Riemann tensor and the covariant derivatives, 
therefore, reads 

L = ai V p R\7 p R + a 2 V r R Pq V r R pq + asVtRpqrsV* R pqrs . (51) 

For general values of a%, a 2 and a 3 , (1511) leads to six order differential equations. We would like to find all values of a%, 
a 2 and a 3 for which (|51[) leads to fourth order equations (imposing the E-criterion on f|5 1 1) is the same as requiring it 
to lead to fourth order equations). Since C 2 (I40|) leads to fourth order equations, and equations of motion are linear 
in term of the Lagrangian density, subtracting a multiplication of £2 from (|5ip does not change the differential degree 
of the equations of motion derived from (fBTj) . We do the following subtraction: 

L = L- a 3 C 2 = (ai - a 3 )V p RV p R + (a 2 + 4a 3 ) V r R pq V r R pq . (52) 

Any values of ai,02 and a 3 that leads to fourth order equations derived from (|5ip. leads to fourth order equations 
derived from (|52|) and vice versa. We define the following constant values in order to write (|52|) in a more compact 
form: 

61 = ai - a 3 , (53) 
b 2 = a 2 + 4a 3 . 

Using the above constant values, L reads 

L = biV p RV p R + b 2 V r R pq V r R pq . (54) 

We will prove that only for b\ = b 2 = 0, L leads to fourth order equations. Then b\ = b 2 = implies that L 2 (up to 
an overall factor) is the most general Lagrangian density, quadratic in term of the Riemann tensor and quadratic in 
the covariant derivative leading to fourth order equations for the metric's components. 

Instead of considering a general metric and calculating the functional variation of (|52[) . let us consider a general 
time-independent spherical metric in d = 4: 

dr 2 

ds 2 = -A(r)dt 2 + -, - , - + r 2 (d9 2 + sin 2 9d(f> 2 ) , (55) 
A{r)B(r) 

and calculate the functional variation of (|54|) with respect to A(r) and B(r). It is known that the functional variation 
does not generally commute with imposing symmetries on the solution; here we have imposed spherical symmetry and 
time translation. We notice, however, our imposed symmetries are isometries of the Riemann manifold -supposedly a 
smooth manifold- so the principle of symmetric criticality is met [24j . In other words the functional variation of the 
action corresponding to (pylj) computed for (|55|) . indeed gives the equations of motion of A(r) and B{r). We will show 
that only for b\ = b 2 = 0, L leads to fourth order equations for A(r) and B(r). So &i = b 2 = must hold in order to 
have fourth order equation for a general metric. 



This term relies within the family of the Lagrangians which are functional of the Riemann tensor not its covariant derivative. In this 
family of the Lagrangian only the ordinary Lovelock gravity meets the E-criterion. 
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In the following we are going to compute the equations of motion for A(r) and B(r) in (|55p. The action corresponding 
to ([54")) computed for (|55|) reads 

S L =4njdt Jdrr 2 ^1(61 V p i?VJ?| Eq + 6 2 V r i? pg Vi?P*| Eq gg) , (56) 
4r 6 \7 a RV a R\ Eqm 



where 



AB 



(-B"A'r 3 - 2,B'A"r z - UB'A'r 2 - 2BA"'r 3 (57) 
- 8BA"r 2 + ABrA' - AAB"r 2 - 8 + 8 AB) 2 , 

— ^ = 9 6 + 8 o r 4 B' A" A' B + 12 r 6 B' A" B A'" - 32 r 3 B 2 A" A' + 8 r 5 B" A' B A" (58) 

AB v ' 

+4 r 6 B" A' B A'" + 64 A B' r + 96 B A" r 2 + 24 r 5 B' A' B A'" - 48 r 3 B' A' 2 B 

+12 r 5 B" A' 2 B' -8r 4 B'AB A"' + 52 r 4 B' A B" A' - 16 r 3 A' B A B" 

-8r 4 B" A' 2 B - 192AB + 12 r 5 B' A" AB" + 8r 5 B A'" A B" 

+32 r 4 B A" A B" - 96 r 2 B' A A' B + 16 r 3 B A" B' A + 36 r 5 B' 2 A" A' 

+6 r 6 B" A' B' A" + 16 r 5 B 2 A'" A" + 24 r 5 B' A" 2 B + 32 A B" r 2 

+112 B' A' r 2 + 96 B 2 A 2 - 32 A 2 B" r 2 B - 64 A 2 B' r B - 96 B 2 A" r 2 A 

+4 r 5 B" 2 A' A - 8 r 3 B' A 2 B" - 12 r 4 B' 2 A A" - 16 r 4 B 2 A'" A' 

+12 r 4 A 2 B" 2 + 20 r 2 B' 2 A 2 + 40 r 4 B 2 A" 2 + 4 r 6 B 2 A'" 2 + 9 r 6 B' 2 A" 2 

The functional variations of (|55)) with respect to A(r) and B(r) at most are functional of the first six derivatives of 
A(r) and B(r), and r. It is straightforward to calculate them. It is then easy to show that 

d f5S~ T \ rB 



9B( 5 )(r 



(ja) = r -^( 2 (^i + b 2 )A + (2b 1+ b 2 )rA'), (59) 



If the equations of motion of A and B are fourth order, then ([5^)1 must algebraically vanish. It algebraically vanishes 
only for 

&i = b 2 = 0. (60) 

Since there exists no non- vanishing value for b\ and 62 leading to fourth order equations for spherical time-independent 
metric, then there exists no non-vanishing value for b\ and b 2 that leads to fourth order equations for components 
of a general metric. To put it in other words, (l60l) besides (|53|) proves that C 2 (|40)) is the only Lagrangian in the 
form of (IBTjl that satisfies the E-criterion. We expect that such a uniqueness property could be generalized to also 
include other C n : C n is the only action constructed from the multiplication of n-times of the first covariant derivative 
of the Riemann tensor which satisfies the E-criterion. In other words, we expect it to be proved that (|3"3"|) is also the 
necessary condition for the E-criterion, and (|36l) is its only solution. 



III. SUMMARY AND OUTLOOKS 



The differential degree for the equations of motion of metric generally is two degrees higher than that of the first 
variation of the Lagrangian with respect to the Riemann tensor if the Riemann tensor were independent of metric (the 
E-tensor). So requiring the same differential degree for the metric's equation and the E-tensor serves as a criterion to 
single out a subset of Lagrangians from a given larger family of Lagrangian. This criterion (the E-criterion) can be 
utilized when one wishes to fix the field redefinition ambiguities [271 ] in quantum loop corrections. 

We have noticed that Lovelock gravity is the result of applying our criterion to Lagrangians which are functional 
of the metric and the Riemann tensor. We have considered Lagrangians in the form of C(gij,^7 a Rijki) polynomial 
in term of the first covariant derivative of the metric. We have found that (|36[) meets the criterion. In particular in 
d = 4, (|4T)]) lead to fourth order differential equation for the metric. So it should be considered when one addresses a 
general fourth order gravity [26j . 
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Perhaps we can apply the E-criterion on other families of Lagrangians: families that includes arbitrary higher 
derivatives of the Riemann tensor. The E-criterion then chooses a subset of these Lagrangians. 

Notice that our criterion is a weaker condition than that of requiring (a strong form of) consistency between Palatini 
and metric formulations [l6l[l7|. It, however, sounds interesting to check if adding terms cubic in the Riemann tensor 
to (|40j) . would lead to a Lagrangian whose Palatini and metric formulations are (strongly) consistent. 
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